On the Formation of Discontinuities of Wave Fronts in a Saturated Granular Body

The paper deals with propagation of plane wave fronts in a solid with a non-linear relation between stress and deformation. The objective is to calculate the distance that a wave front covers before it loses continuity. The formulae derived for a general quasi-linear system of two equations are applied to the propagation of plane compression waves in dry and partially saturated granular bodies. In the case of a saturated body with gas bubbles, the influence of gas and capillary pressure on the stiffness of the body is taken into account. Numerical calculations relevant to soil mechanics are presented. For the numerical calculations a constitutive equation of the hypoplasticity theory for granular materials has been used. Received November 1, 1997     

Dynamic analysis of geometrically nonlinear robot manipulators

In this paper, a method for the dynamic analysis of geometrically nonlinear elastic robot manipulators is presented. Robot arm elasticity is introduced using a finite element method which allows for the gross arm rotations. A shape function which accounts for the combined effects of rotary inertia and shear deformation is employed to describe the arm deformation relative to a selected component reference. Geometric elastic nonlinearities are introduced into the formulation by retaining the quadratic terms in the strain-displacement relationships. This has lead to a new stiffness matrix that depends on the rotary inertia and shear deformation and which has to be iteratively updated during the dynamic simulation. Mechanical joints are introduced into the formulation using a set of nonlinear algebraic constraint equations. A set of independent coordinates is identified over each subinterval and is employed to define the system state equations. In order to exemplify the analysis, a two-armed robot manipulator is solved. In this example, the effect of introducing geometric elastic nonlinearities and inertia nonlinearities on the robot arm kinematics, deformations, joint reaction forces and end-effector trajectory are investigated.     

Theoretical and experimental analysis of longitudinal wave propagation in cylindrical viscoelastic rods

Wave propagation in viscoelastic rods is encountered in many applications including studies of impact and fracture under high strain rates and characterization of the dynamic behavior of viscoelastic materials. For viscoelastic materials, both material and geometric dispersion are possible when the diameter of the rod is of the same order as the wavelength. In this work, we simplify the Pochhammer frequency equation for low and intermediate loss viscoelastic materials and formulate corrections for geometric dispersion for both the phase velocity and attenuation. The formulation is then experimentally verified with measurements of the phase velocity and attenuation in commercial polymethylmethacrylate rods that are 12 and in diameter. Without correcting for geometric dispersion, the usable frequency range for determining the phase velocity and attenuation for the rod is about , and about for the rod. Using the correction procedure developed here, it was possible to accurately determine the phase velocity and attenuation up to frequencies exceeding for the rod and for the rod. These corrections are applicable to many polymers and other viscoelastic materials. From thereon, the viscoelastic properties of the material can be determined over a wide range of frequencies.     

Regularity and asymptotic behavior of solutions of nonautonomous differential equations

The existence of a nonautonomous approximate inertial manifold is shown for problems of the formu + Au + N(t,u)=0, in whichA is a self-adjoint operator with compact resolvent in a Hilbert spaceH. The operatorN(t, u) = G(u) + F(t, u) is nonlinear withG a monotone gradient that is locally Lipschitz fromD(A ^1/2) intoH, andF:⁺×HH a Lipschitz perturbation that is Hölder continuous int. Weak solutions are shown to be uniformly locally Hölder continuous intoD(A) with equicontinuity in families of solutions with ¦u(0)¦ r.A priori estimates of ¦Au(t)¦ are also verified and used in a skew-product flow to show there is a global attractor whose component elements form a equicontinuous family of solutions.     

Dimensional scaling of flame propagation in discrete particulate clouds

The critical dimension necessary for a flame to propagate in suspensions of fuel particles in oxidiser is studied analytically and numerically. Two types of models are considered: First, a continuum model, wherein the individual particulate sources are not resolved and the heat release is assumed spatially uniform, is solved via conventional finite difference techniques. Second, a discrete source model, wherein the heat diffusion from individual sources is modelled via superposition of the Green's function of each source, is employed to examine the influence of the random, discrete nature of the media. Heat transfer to cold, isothermal walls and to a layer of inert gas surrounding the reactive medium are considered as the loss mechanisms. Both cylindrical and rectangular (slab) geometries of the reactive medium are considered, and the flame speed is measured as a function of the diameter and thickness of the domains, respectively. In the continuum model with inert gas confinement, a universal scaling of critical diameter to critical thickness near 2:1 is found. In the discrete source model, as the time scale of heat release of the sources is made small compared to the interparticle diffusion time, the geometric scaling between cylinders and slabs exhibits values greater than 2:1. The ability of the flame in the discrete regime to propagate in thinner slabs than predicted by continuum scaling is attributed to the flame being able to exploit local fluctuations in concentration across the slab to sustain propagation. As the heat release time of the sources is increased, the discrete source model reverts back to results consistent with the continuum model. Implications of these results for experiments are discussed.     

Various solution methods,accompanied by dynamic investigation,for the same problem as a means for enriching the mathematical toolbox

A geometrical task is presented with multiple solutions using different methods, in order to show the connection between various branches of mathematics and to highlight the importance of providing the students with an extensive ‘mathematical toolbox’. Investigation of the property that appears in the task was carried out using a computerized tool.     

Geometric phase of an open double-quantum-dot system detected by a quantum point contact

We study theoretically the geometric phase of a double-quantum-dot(DQD) system measured by a quantum point contact(QPC) in the pure dephasing and dissipative environments, respectively. The results show that in these two environments, the coupling strength between the quantum dots has an enhanced impact on the geometric phase during a quasiperiod. This is due to the fact that the expansion of the width of the tunneling channel connecting the two quantum dots accelerates the oscillations of the electron between the quantum dots and makes the length of the evolution path longer.In addition, there is a notable near-zero region in the geometric phase because the stronger coupling between the system and the QPC freezes the electron in one quantum dot and the solid angle enclosed by the evolution path is approximately zero,which is associated with the quantum Zeno effect. For the pure dephasing environment, the geometric phase is suppressed as the dephasing rate increases which is caused only by the phase damping of the system. In the dissipative environment,the geometric phase is reduced with the increase of the relaxation rate which results from both the energy dissipation and phase damping of the system. Our results are helpful for using the geometric phase to construct the fault-tolerant quantum devices based on quantum dot systems in quantum information.     

Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.     

Failure Impact,Availability and Cost Analysis of PONs Based on a Network Geometric Model

Passive Optical Networks(PONs)are considered as the preferred solution for broadband fibre-based access networks.This is because PONs present low cost deployment,low energy consumption and also meet high bandwidth demands from end users.In addition,end users expect a high availability for access networks,while operators are more concerned about reducing the failure impact(number of clients affected by failures).Moreover,operators are also interested in reducing the cost of the access network.This paper provides a deep insight into the consequences that the physical topology and design decisions cause on the availability,the failure impact and the cost of a PON.In order to do that,the physical layout of the PON deployment area is approximated by a network geometric model.A PON deployed according to the geometric model is then assessed in terms of failure impact,availability and cost.This way,the effects of different design decisions and the physical layout on these three parameters are evaluated.In addition,the tradeoffs between availability,failure impact and cost caused by planning decisions and the physical topology are identified and pinpointed.     

Thermodynamic description of the viscosity of liquid solder alloys with minor Co impurities

Because of the increasing complexity and cost of experiments carried out, the data for the multi-component alloy systems have frequently been obtained by numerical modelling. It is clear that the related calculations require reliable data dealing with the pure components and binary alloy systems. Selecting the reliable data concerning the pure components from the literature, the viscosities for the SAC and (SAC)_1?x Co_x solder alloys have been calculated using different viscosity models (geometric and physical). The viscosity decreases as the amount of tin content increases in the SAC387 alloy while the addition of the cobalt to SAC387 solder results in the increasing of the viscosity. Moreover, by computing the root mean square values between theoretical and experimental viscosities, it can be concluded that the lowest value among all models is that of obtained by Kaptay equation.